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ABSTRACT 

We examine the optical emission line properties of Brightest Cluster Galaxies 
(BCGs) selected from two large, homogeneous datasets. The first is the X-ray selected 
National Optical Astronomy Observatory Fundamental Plane Survey (NFPS), and 
the second is the C4 catalogue of optically selected clusters built from the Sloan 
Digital Sky Survey Data Release 3 (SDSS DR3). Our goal is to better understand 
the optical line emission in BCGs with respect to properties of the galaxy and the 
host cluster. Throughout the analysis we compare the line emission of the BCGs to 
that of a control sample made of the other bright galaxies near the cluster centre. 
Overall, both the NFPS and SDSS show a modest fraction of BCGs with emission 
lines (~ 15%). No trend in the fraction of emitting BCGs as a function of galaxy mass 
or cluster velocity dispersion is found. However we find that, for those BCGs found 
in cooling flow clusters, 71^ 4 have optical emission. Furthermore, if we consider only 
BCGs within 50kpc of the X-ray centre of a cooling flow cluster, the emission-line 
fraction rises further to 100^5 %. Excluding the cooling flow clusters, only ~ 10 % of 
BCGs are line emitting, comparable to the control sample of galaxies. We show that 
the physical origin of the emission line activity varies: in some cases it has LINER- 
like line ratios, whereas in others it is a composite of star-formation and LINER-likc 
activity. We conclude that the presence of emission lines in BCGs is directly related 
to the cooling of X-ray gas at the cluster centre. 

Key words: galaxies: clusters, cooling flows:general - galaxies: evolution - stars: 
formation - galaxies: stellar content - surveys 



1 INTRODUCTION 

The brightest cluster galaxy (BCG) is typically a giant, 
red elliptical or cD galaxy, located near the centre of 
the gravitational potential. It is likely that a rich history 
of galaxy-galaxy interactions and mergers is responsible 
for the unique morphology of such galaxies. This is sup- 
ported indirectly by several pieces of evidence, including: 
the lu minosity of the BCG is corre lated with the cluster 
mass (IScottJll957l: iLin fc Mohrll2004h and X-ray luminosity 
{Hudson fc Ebelingl 11997^ 7 BCGs in the most X-ray lumi- 
nous clusters are larger and have surface brightness profiles 
which ar e less steep than the ir low X-ray luminosity coun- 
terparts (|Brough et al l l2005l h and the velocity dispersion 
of BCGs rises l ess steeply with luminosity than for other 
bright galaxies (|von der Linden et al.| [2006h All this indi- 
rect evidence for a rich merger history is supported by high 



resolution imaging of BCGs obtained with the Hubble Space 
Telescope, which has revealed that the cores of these galaxies 
can be complex, ofte n showing multip le nuclei and promi- 
nent dust signatures (|Laine et al.ll2003t ). 

However, the evolutionary history of these galaxies is 
still not completely understood, and work on cooling flow 
(CF) clusters has hinted at another possible mechanism 
for adding to the stellar mass of BCG. The original cool- 
ing flow hypothesis is that hot cluster X-ray gas cools 
and condenses out of the intracluster medium (ICM) into 
the cl uster's poten tial well, forming molecular clouds and 
stars (Fabian 1994). This drop-out occurs at the centre of 
the cooling flow, within the cooling radius, ie. onto the 
BCG. Cooling flow clusters are common in the local uni- 
verse (making up abo ut 50% of the population in an X - 
ray flux limited sample IPeres et al.|[i"99sl ; IChen et alj|2007l )'). 
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and cD galaxies are often found at the centre of these 
systems. Because of this, a link between the cooling X- 
ray gas and recent star fo rmation in th e BCG has been 
discussed for many years (|Fabian|[l993) . Convincing ob- 
servations which support thi s idea have been presen ted: a 
blue and UV-colour exce ss ijMcNamara et al.1 1 19961. 
Hicks fc Mushotzkv 20051). molecular gas (|Edge et al.l 
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Jaffc et al. 2001; Salome & Combes 2003) and Ha emission 



ilCrawford et al.lll999l ; iDonahue et al.ll200d ; ICrawford et al.1 
2005) have all been seen in the cooling flow cluster's 
BCG. However, although the morphology of the Ha emis- 
sion is diffuse and filamentary, indicating star formation 
in some CF BCGs, in others it is very co mpact and more 
chara c teristic of AGN dominated emission (IDonahue et all 
200d; lEdwards fc Robert! 120071 : iHatchl 120071 ). As well. 



von der Linden et al.l 1 200(: ) recently showed that optical 



emission lines in BCGs predominantly arise from LINER 
emission, rather than normal star formation. 

More recent X-ray satellite measurements from Chan- 
dra and XMM have shown that the gas does not cool di- 
rectly from th e hot X-ray phase thr ough to the cool molecu- 
lar gas phase l|Bohringer et al.ll2002T ). but rather that a large 
amount of the cooling gas is being reheated before condens- 
ing out of the ICM. The current paradigm is that AGN ac- 
tivity in the BCG is reheating the cooling X-ray gas, which 
implies a more complicated feedback pr ocess between the 
coolin g gas and the central galaxy (e.g Pizzo lato fc Sokeij 
2005). This leads to revised, predicted mass deposition rates 
that are now in reasonable agreement with the obs erved val- 
ues (|Bohringer et al.ll2002l ; |pizzolato fc Sokerl2005l ). The ob- 
served molecular and ionic gas may be attributed to a small 
amount that has cooled from the cooling flow; alternatively, 
the Ha may be excited by the AGN itself. Detailed studies 
of star formation indicators in galaxy groups and clusters in 
cooling flows and non-cooling flows, discriminating between 
those with and without AGN activity, are required in order 
to analyze the relative importance of the different processes. 

A correlation between optical line emission in the BCG 
and cluster pr operties has bee n explo red by several authors, 
most notably ICrawford et al.l (|l999l ). who found that 27% 
of BCGs have optical line emission, and that the projected 
distance from the BCG to the X-ray centre is less for line 
emitting galaxies than for n on-e mitting galaxies. Recently, 
Ivon der Linden et al.1 (|2006l ) and lBest etaLI (|2006l ) have ex- 
plored the propert ies of BCGs in th e SDSS, using the C4 
cluster catalogue (Mil ler et al 1 120051 ). These authors find 
that radio-loud AGN activity is more frequent in BCGs than 
in other galaxies of the same mass, but that this frequency 
does not depend strongly on cluster velocity dispersion. On 
the other hand, in their study of radio-loud p r opert ies of 
an X-ray selected cluster sample, iLin fc Mohrl (2006) find 
the overall radio-loud fraction to be 30% in BCGs, and that 
the fraction is higher i n mor e massive clusters. Importantly, 
Ivon der Linden et al.l (|2006l ) find that many of these radio- 
loud galaxies would not necessarily be identified as AGN 
from their optical emission lines and, in fact, that optical 
AGN activity appears to be less frequent among BCGs than 
other cluster galaxies of similar mass. The interpretation is 
complicated by the fact that the radio-selected galaxy sam- 
ple, though restricted to red galaxies, could be contaminated 
by galaxies in which the low luminosity radio emission arises 
from star formation, rather than AGN activity. 



In this paper, we explore optical line emission in BCGs 
with respect to properties of the galaxy and the host clus- 
ter, using two large, homogeneous datasets. One sample is 
taken from the X-ray selected, National Optical Astronomy 
Observatory Fundamental Plane Survey (NFPS), for which 
the X-ray properties are known for all clusters. For many 
of these clusters, we are able to identify those with short 
cooling times (CF clusters) based on ROSAT, Chandra, 
or XMM-Newton observations. We complement this sam- 
ple with optically-selected clusters drawn from the Sloan 
Digital Sky Survey Data Release 3 (SDSS DR3), which is 
not biased toward X-ray luminous clusters, and is therefore 
more representative of the cluster population. In addition, 
the greater spectral coverage of the SDSS allows us to use 
emission line ratios to identify whether the emission arises 
predominantly from composite HII region and LINER ac- 
tivity, or from LINER activity alone (pure HH-region, and 
Seyfert-like emission are both rare). 

The paper is organized as follows. In section [2j we in- 
troduce our galaxy samples and selection criteria. In section 
[3] we report our results. For the X-ray selected NFPS we 
compare the frequency of H/3 emission in BCGs as a func- 
tion of BCG magnitude and distance to the cluster centre. 
Similar results are found for the SDSS sample, based on the 
Ha emission line. However, we find that there are differ- 
ences between the BCGs in the two samples, which we can 
attribute to the nature of the X-ray emitting gas. In section 
[4] we consider the impact of our results on various galaxy 
and cluster formation hypotheses, and summarize our re- 
sults. We conclude in section [S] Unless otherwise stated our 
analysis assumes the values f2 m =0.3, for the matter den- 
sity parameter, £7a=0.7 for the cosmological constant, and 
-Ho=100 km/s/Mpc for the Hubble parameter. Lx refers to 
the bolometric X-ray luminosity throughout. 



2 DATA AND SAMPLE SELECTION 

For both NFPS and SDSS, our goal is to identify the BCG 
and a similar sample of luminous "control" galaxies that is 
located in the inner regions of the cluster. We first discuss 
the sample selection of the NFPS. 



2.1 NOAO Fundamental Plane Survey 

The NFPS is an all sky study of 93 X-ray select ed rich clus- 
ters with redshifts between 0.010 and 0.067 l|Smith et al.l 
l2004h . The goals of the project are to measure cosmic flows 
on the scales of 100 h _1 Mpc and to build a large homo- 
geneous sample with which to investigate physical effects 
and environmental i nfluences on early - type galaxy evolu- 
tion and formation l|Smith et al.l [2004) . The spectroscopic 
observations are made though a fiber diameter of 2" and are 
limited to red sequence galaxies; galaxies more than 0.2 mag- 
nitudes bluer than the red sequence were not generally ob- 
served spectroscopically. There are spectra for 5388 galaxies 
with a wavelength cov erage between 400 and 6100 A and 
a resolution of ~ 3 A l|Smith et al.ll2004l ). 
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2.1.1 Cluster, BCG, and Control Sample Definitions 

In order to respect the completeness of the NFPS clus- 
ter sample, we exclude the 13 clusters that were observed 
serendipitously and that did not meet the original Lx limit 
of 10 42 ' 6 erg/s. For each cluster, the redshifts were used to 
calculate the cluster velocity dispersion and the radius at 
which the mean density interior to the cluster is 200 times 
the critical density (er c ; and r2oo, respectively). We used 
the prescription r2oo= y/So' c i /WOO km/s/Mpc as derived in 
ICarlberg et ail l|l997t) . Note that r2oo is typically of order 
1.5 h~ x Mpc, much larger than the cooling radius of a typ- 
ical cooling- flow cluster (about 200 kpc). The centre of the 
cluster is taken to be at the peak of the X-ray emission 
labeling et al.lll996l . |2000| ). The galaxies are then assigned 
to th e clusters ba s ed on a radial and a velocity weighting 
as in ISmith et al.l ([2004). For this analysis, we are inter- 
ested in only the bright galaxies in the central regions of 
the cluster and so include a magnitude limit of Mk < —24, 
based on K-band total magnitu des obtained from 2MASS 
catalogues (|Skrutskie et al.ll2006l ). This is about half a mag- 
nitude brighter t han the charact eristic K-band magnitude of 
M K * = -23.55 (|Lin et al.ll200j . for H = 100 km/s/Mpc). 
Since we are interested only in galaxies occupying a similar 
environment as the BCG, ie. in the central regions of the 
cluster, we consider only galaxies within 0.5r2oo, and with 
velocity differences with respect to the cluster mean velocity 
less than twice the cluster velocity dispersion. 

The BCG is then defined as the first rank cluster galaxy 
using the K-band magnitudes. In twenty cases, the BCG was 
not observed spectroscopically due to constraints in the fiber 
positioning, and these clusters have been excluded from our 
analysis. Our final sample consists of 60 clusters with BCGs, 
as summarized in Table Q] and listed in Tabled Our control 
sample consists of the 159 other bright (Mk < —24) galaxies 
within the same radial and velocity cuts described above. 
The BCGs are, of course, excluded from the control sample. 

For the NFPS, we use the stellar absorption-corrected 
H/3 emission as an indicator of star formation activity, since 
Ha is not generally available in these spectra. Although Hp 
emission is relatively weak, the high signal-to-noise ratio of 
these spectra allow us to measure its strength reliably af- 
ter correcting for underlying stellar absorption (see the er- 
rors quoted in Table [2}. The observed galaxy spectra are 
divided by best fit a bsorption templ ate stellar population 
synthesis models from lVazdekisI (|l999h . which have been red- 
shifted and broadened to match the velocity dispersion of the 
observed galaxy. Subsequently, the H/3 equivalent width is 
measured directly, without assuming a particular line pro- 
file, from the ratio of the o bserved spectrum to the best 
fit model l|Nelan et al.ll2005l ). A thorough disc ussion of the 
emissi on line measurements can be found in iNelan et al.l 
|2005l ). For two of the NFPS BCGs, in Abell 780 and in 
Abell 1795, the nebular emission is strong enough that our 
standard methods for obtaining reliable velocity dispersions, 
and hence stellar absorption corrections, fail. For galaxies of 
similar magnitude, the typical stellar absorption at H/3 is 
~1.5 A with an uncertainty of 0.3 A. This error dominates 
our uncertainity in the total equivalent width of emission 
for these two special cases. For the BCGs in Abell 780 and 
Abell 1795, the H/3 equivalent widths are 7.8 A and 7.2 A, 
respectively. However, to achieve a plot that is more eas- 



ily read, these points are set to respective lower limits of 
3.3 A and 3.7 A in Fig. [5] and Fig. H 

We will define emission-line galaxies to be those with an 
equivalent width > 0.5 A. Since the Ha and [Nil] lines are 
genera lly unavailable, we are unable to use iBaldwin et alj 
l|l98ll . hereafter BPT) dia grams to reliably distinguish emis- 
sion due to star formation from that arising from AGN ac- 
tivity. 



2.1.2 Cooling Flow definition 

In general, we designate an NFPS cluster as "cooling flow" 
(CF) or "non-cooling flow" using mass depositi on rates and 
X-ray cooling time s from published ca t alogues (IPeres et alj 
1 19981 ; IWhitel |2000| ; I Allen et al.l l200ll ; iBirzan et al.l |2004| ). 
This is a somewhat subjective classification complicated 
by the fact that not only are mass deposition rates cal- 
culated from ROSAT observations typically 2-10 times 
higher than those cal culated from Chandra observations 
jBohringer et alj|2002l ), but also higher resolution spectra 
from XMM Ne wton are not well m atched to an isobaric 
cooling model l|Peterson et alj [20031 ); thus the mass depo- 
sition rate may not be an exact indicator of a cooling flow 
cluster. Therefore, whenever possible, we prefer to use re- 
cent cooling flow designations based on the presence of a 
central temperature gradient in XMM Newton or Chandra 
observations. For the rest, we are left with using the mass 
deposition rate based on ROSAT data as an indicator. For 
the 11 cases where observations are from the Chandra or 
XMM Newton satellites, we define a CF cluster to be one 
with a mass deposition rate M > 0; otherwise, we require 
M > lOOMQyr -1 . Within this framework we have 14 CF 
clusters and 19 non-CF clusters. The CF status of another 
27 clusters is unknown. Clearly this is not an unassailable 
definition; however, it is likely that most, if not all, of the 
clusters we classify as CF clusters really do have short cool- 
ing times in the centre. Some galaxies with low mass de- 
position rates will undoubtedly fall in our non-CF sample. 
With this in mind, our results are not sensitive to this def- 
inition. We further discuss this point in Section \3. II as well 
as a continuous method of defining a CF cluster, based on 
an excess of observed X-ray luminosity to predicted values 
([McCarthy et alj|2004l ). 

2.2 Data from the Sloan Digital Sky Survey 

Our second sa mple of galaxy cl usters is derived from the 
C4 catalogue (|Miller et al.ll2005h . based on the third data 
release (DR3) of the SDSS. This release covers 5282 square 
degrees in imaging, and 4188 square degrees in spectroscopy 
l|Abazaiian et al1l2005T ). Imaging was taken in five optical 
bands, u', g', r', i' and z' with a median spatial resolution 
of 1.4" in r'. Spectra are observed with a aperture diam- 
eter of 3", cover the wavelength range from 3800 to 9200 
A and have a spectral resolution of ~3 A in r'. Redshifts 
measured from these s pectra are accurate to ~ 30 km/s 
l|Abazaiian et all 120051 ). The Ha line strengths are mea- 
sured by fitting gaussians to th e line profile in a standard 
pipeline ([Stoughton et all 120021 ). For W (Ha) > 5A, the 
equivalent width uncertainty is less than 20%. For weaker 
lines, Wo(Ha) < 5A the errors are known to be large 
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l|G6mez et alj|2003h . We have made no correction to the 
emission-line strengths for underlying stellar absorption. For 
our purposes this is safe to neglect because, whereas even 
a modest star formation rate generates considerable Ha 
emission, the stellar absorption doesn't vary by more than 
~ 1 A for moderate ly old populations. 

As described in iMiller et all (|2005l ) , clusters and groups 
are identified as overdensities in the multi-dimensional space 
spanned by position, redshift, and five-colour photometry. 
There are 1106 clusters identified in the SDSS DR3 using 
this algorithm. We further select objects with redshifts z < 
0.10, to minimize the incompleteness of the sample. In or- 
der to reduce the number of clusters with uncertain velocity 
dispersions, we exclude clusters flagged as having significant 
substructure. Specifically we include only clusters flagged as 
SUB=0, ie. those for which the ratio of the standard devi- 
ation of the cluster velocity dispersion profile to the mea n 
cluster velocity dispersion is less than 15 l|Miller et alj2005l ). 
This reduces the number of clusters in our sample to 825. 



2.2.1 Cluster, BCG, and Control Definitions 

For the SDSS we define the BCG and clusters in as close a 
manner as possible to the NFPS case. Again, the BCG is 
the brightest galaxy in the K band (and with Mk < — 24), 
within half of r2oo of the geometric centre of the cluster, 
and within two times the cluster velocity dispersion. How- 
ever, the geometric centres we use are different than those 
in the C4 catalogue. We start with the C4 catalogue geo- 
metric centres measured using the luminosity weighted po- 
sition average of all the galaxies within IMpc and four times 
the velocity dispersion. As this definition encompasses such 
a large area, multiple substructures along the line of sight 
heavily influence the position of the centre. For example, 
the geometric centre will be placed in between two obvious 
subclumps. Thus, we iteratively recalculate a luminosity- 
weighted centre using only galaxies within two times the 
cluster velocity dispersion, and with magnitude Mk < —24. 

Due to the incomplete spectral coverage of the central 
galaxies, there are 154 cases in which the brightest cluster 
galaxy had not been observed spectroscopically. We exclude 
these clusters from the sample. To ensure this is not in- 
troducing an important bias, we verified that the subset of 
BCGs observed in the photometric catalogue but not in the 
spectroscopic catalogue have an equivalent (it' — r') colour 
distribution to those in our final, spectroscopic sample. Fi- 
nally, we also remove those clusters whose geometric centres 
are within 15' of a survey boundary, as well as those with 
measured velocity dispersions greater than 1200 km/s. The 
latter restriction is made because such high values are usu- 
ally a result of significant contamination from line-of-sight 
substructure. These selections leave us with a final sample 
of 328 BCGs. 

Our control sample is built from the other bright galax- 
ies near the centre of the clusters, as in the NFPS case. 
In order to increase the size of our control sample, we in- 
clude clusters where the brightest central galaxy was not 
measured spectroscopically. There are 526 control galaxies 
in 353 groups and clusters with a velocity dispersion less 
than 1200km/s. We note that the number of control galax- 
ies per cluster is significantly larger in the X-ray luminosity 



selected NFPS clusters, probably because they are richer 
than the optically-selected SDSS clusters. 

Because the centres of the NFPS clusters are based 
on the X-ray centroid, whereas in the SDSS the geomet- 
ric centre is used, it is important to compare the two 
definitions. In order to find X-ray centers for the SDSS 
clusters, we have match ed them to X-ray cl uster cata- 
logues including: NOR AS feohringer et al. 2000), REFLEX 
l|Bohringer et aLlbOOlT) , BAX dSadat et al.1 120041) XBA CS 
lEbeling et al.lll996ri. RASS C ALS llMahdavi et alll2000D. as 



well as 



Popesso et al] (I2004T ); iMulchaev fc Zabludofj (|l998h 



and lHorner et al.l I 200lf ). Most of the cluster catalogues are 
based on ROSAT observations, where the flux limit is high. 
This restricts our sample of SDSS clusters with X-ray de- 
tections to the small subset of massive clusters at z < 0.03. 
There are 35 X-ray cluster matches, ie. cases where the X- 
ray centre is within half of r2oo of one of our 328 SDSS 
clusters. Fig. [T] shows that for most of the matched clus- 
ters, there is good agreement between the X-ray and geo- 
metric centres. However, also depicted in the figure is how 
some cases exhibit differences of up to 0.5 Mpc, especially for 
lower mass clusters. It is unlikely that this is caused by an 
uncertainty in the X-ray centres, as the centres of the NO- 
RAS, BCS, and XBACS samples were found by determining 
the two-dimensional centre of mas s using the Voronoi Tes- 



sellat ion and Percolation method (|Ebeling fc Wiedenmann 
1993), and are generally accurate to about 1' ( Ebeling et al. 
Il998f ). corresponding to ~ 90 kpc at z ~ 0.08. More likely, is 
that the geometric centers do not trace the gravitation po- 
tential of the cluster as well as the X-ray centres. We note 
that for many of the cases in which the centres are discrepant 
by > 200 kpc, the geometric centre appears to be contam- 
inated by in-falling groups; which is not surprising as it is 
measured using a typical line of sight projection of ~ 3 Mpc 
(for a = 500km/s). 



2.2.2 Emission Line diagnostics 

Optical emission line galaxies in the SDSS are identified from 
the Ha line, which is the line in our wavelength range that 
is most sensitive to star formation activity. To ensure a fair 
comparison with the NFPS, we must choose a threshold in 
Ha equivalent width that is comparable to the H/3 limit used 
in that survey. Recall that the lines are measured using dif- 
ferent techniques, from spectra of different resolution and 
signal-to-noise, and obtained with different fibres; further- 
more, the NFPS H/3 measurements are corrected for under- 
lying stellar absorption, while the SDSS Ha lines are not. 
Therefore we opt for an empirical "calibration" between the 
two, by plotting the NFPS H/3 equivalent widths against the 
SDSS Ha equivalent widths, for galaxies that appear in both 
surveys (Fig. [SJ. Although there is plenty of scatter, there 
is a strong correlation between the two lines, and we find 
the Ha index is about four times larger than the H/3; en- 
couragingly, this is comparable to the factor of 4.5 derived 
by comparing the Ha/H/3 ratio for the subset galaxies in 
the NFPS for which measureme nts of both emission lines 
are available IjNelan et alj I2005T ) . Thus, an Ha equivalent 
width cut of 2 A is comparable to an H/3 equivalent width 
of 0.5 A, and the fraction of galaxies in our samples above 
either of these thresholds is similar. Note that the points 
at Wo(Ha) — are non-detections (plotted with the aver- 
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Figure 1. Cluster velocity dispersion as a function of the differ- 
ence in the X-ray and geometric centres for the 35 SDSS BCGs 
that have available X-ray positions. The clusters with large dis- 
crepancies between the two centres are generally found to have 
geometric centres highly influenced by in-falling groups. 

age uncertainty of 0.15A), and those with W (Ha < 0) are 
detected in absorption. We have experimented with SDSS 
Ha equivalent width cuts of 2±1 A and find that our fi- 
nal results do not change significantly. The good correlation 
also gives us additional confidence in the template correction 
used to correct the H/3 equivalent widths for stellar absorp- 
tion, which is relatively much more important here than for 
Ha (the average stellar absorption strength is ~ 1 A, with 
variations at the ~ 0.15 A level). 

Due to the proximity of the galaxies in our sample and 
the finite fiber size though which they are observed, the 
amount of extended line emission could be underestimated. 
However, as we do not know if the emission is extended we 
have not explicitly accounted for the finite fiber size, nor for 
the difference in fiber diameters between two surveys. Here 
again, we argue that the effect on the line strength measure- 
ments is calibrated by our use of an empirical relation. 

Ha line emission may arise from either ionization by 
hot stars, or ionization from AGN activity. We use the 
AGN classifi cation taken from the emission-line analysis 
discussed in iKauffmann et all (|2003T) . which is based on 
Ha/ [Nil] and [Om]/H/3 diagnostic ratios, which separates 
AGN, star- forming, and composite (intermediate between 
star-forming and AGN) regions on BPT diagrams. Table 
[1] shows that, in our SDSS sample, ~ 65 % of galaxies 
with emission have line ratios consistent with AGN or com- 
posite emission, and this fraction may be somewhat higher 
for th e BCG p o pulat ion, relative to the control galaxies. 
Using IHo et all (Il997f> to separate Sey ferts from LINERs, 
along with the Kauffmann et al. (2003) definitions, we find 
that ~ 33% of emitting SDSS BCGs can be reliably mea- 
sured as LINERs, ~ 27% as composite; Seyfert-like emission 
is negligible ( ~ 3%). If we assume the ionizing source for 
the HII regions is stellar, and that all of the Ha emission 
is within the fibre diameter, then typical Ha luminosities 
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Figure 2. The absorption-corrected NFPS H/3 equivalent width 
versus the SDSS Ha equivalent width (uncorrected for absorp- 
tion) for galaxies present in both surveys. The points with 
W (Ha) = are non-detections (plotted with the average Ha 
error of 0.15 A), and those with W D (Ha) < are detected in ab- 
sorption. We identify line-emitting NFPS galaxies as those with 
Wo (H/3) > 0.5 A. Correspondingly, we use a value of W (Ha) 
> 2 A for the SDSS galaxies. The best fit lines are constrained 
to go through (0,0) and are only fit to galaxies with Wo (Ha) > 
2 A and ignore an outlier at W (Ha) f» 40. The two different 
lines are fitting for Wo (H/3) as a function of Wo (Ha) or vice 
versa. 



of our emitting BCGs co rrespond to star formation rates of 
~ 0.5 to ~ 1.6 Mgyr" 1 i|Kennicutdll99Sh . 



2.2.3 Galaxy Colour 

Fig. [3] shows the distribution of Ha eq uivalent width as a 
functi on of (it' — r') colour. Following IStrateva fc Knappl 
i|200ll ). we use the value of 2.2 for the (u' — r') colour 
separation of blue and red galaxies. As the NFPS galaxies 
are red selected, we include a corresponding colour cut of 
{v! — r') > 2.2 for the SDSS galaxies when directly compar- 
ing to the NFPS results (as in Section [3} . This colour cut 
excludes 10 of the 328 SDSS BCGs (only one with X-ray 
data) as well as 53 control galaxies. 

As seen in Fig. |3j most BCGs are redder than this cut 
whether or not they show optical emission. Thus, if the line 
emission is due to star formation, it does not dominate the 
global colours of these BCGs, which are presumably quite 
old. Because of this, results based on BCGs in the NFPS are 
not likely to be affected by the colour cut in that sample. 
For the control sample, most of the emission-line galaxies 
are bluer than our cut. Therefore we expect the fraction of 
control galaxies with emission, as presented in Section [3] 
to be biased low relative to the BCGs. This is evident from 
TableQ] in the red-selected SDSS sample, the overall fraction 
of emission line galaxies is similar (~ 11 %) for both control 
galaxies and BCGs, while for the unrestricted sample the 
fraction of control galaxies with emission is ~ 18%. This 
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Figure 3. The logarithm of Wo (Ha) as a function of (u' — r') 
colour of the galaxies in the SDSS. The crosses are the control 
sample, and the filled circles are the BCGs. The colour break of 2.2 
is shown as the vertical dotted line, and the horizontal dashed line 
separates the emitting galaxies from the quiescent. The galaxies 
at Log[W (Ha)] = -1.7 are those with Ha absorption. 



Figure 4. The cluster velocity dispersion as a function of the 
bolo metric X-ray l umin osity. The solid line is the preheated model 
from lBabul et al.l [I2002T), which is known to provide a good matc h 
to non-CF clusters jMcCarthv et al.l |20o4 iBaloeh et alll200rj) . 
Known CF clusters in our sample do generally lie to the right of 
this model. 



difference is not large enough to affect any of our conclusions 
in Section O however. 



3 RESULTS 

In this section we examine the line emission in the galaxies 
in our cluster samples as a function of X-ray luminosity, K- 
band magnitude, and distance from the centre of the cluster. 
The numbers of emission-line galaxies in both surveys are 
given in column 4 of Table [T] There is a higher fraction of 
emitting galaxies among the NFPS BCGs; 20±6% are emit- 
ting, compared to only 9±2% of the controls. If we look only 
at those BCGs identified with a CF cluster we find an even 
higher fraction of them show emission, 7li? 4 , while the con- 
trol sample is unchanged (14±6%). On the other hand, only 
11±2% of the red BCGs selected from the SDSS sample 
show emission, comparable to that of the control popula- 
tion (11±2%). In this section we will explore these trends, 
and differences between the two samples, in more detail. We 
use errors derived from the posterior probability distribution 
where sample sizes are small. 



3.1 Dependence on the Presence of a Cooling 
Flow 

The most prominent result from the NFPS clusters is that 
the presence of a cooling flow is highly correlated with the 
presence of emission lines in the BCG. Fig. [4] shows the 
bolometric X-ray luminosity against the cluster velocity dis- 
persion (a proxy for dynamical mass). We can see that most 
of the cooling flow clusters have larger X-ray luminosities 
for their mass, and that it is these clusters that have a BCG 



with line emission. Therefore, in the rest of the NFPS anal- 
ysis we separate our sample into CF and non-CF subsets. 
Notice that in Fig. [S] and Fig. [5] (which we discuss below) , 
none of the emitting BCGs are non-CFs. As mentioned in 
Section \2. 1.1 1 we use a rather strict cut for our definition of 
a CF, requiring ROSAT based M > lOOMQj/r" 1 . However, 
changing our arbitrary definition of a CF cluster (e.g. to 
those with M > lOA^Qyr -1 ) does not significantly change 
our results. 

It would be useful to have a way to identify likely 
cooling-flow clusters without the need for high qual- 
ity surface brightness and temperature maps. Recently, 
iMcCarthv et al.l (|2004h has shown that CF clusters show sig- 
nificantly higher total X-ray luminosities, relative to their to- 
tal dynamical mass, consistent with predictions from steady- 
state cooling models. On the other hand, non-CF clusters 
can be well modelled with cosmological haloes in which 
the gas has been preheated to ~ 300 — 500 keV cm 2 
i|Balogh et al.ll2006| . see also). This suggests that one could 
use the excess X-ray luminosity relative to these preheated 
models as an indicator of CF status. The solid line in 
Fig. 3] repre s ents m odels for clusters with preheating from 
iBabul et alJ < |2002h . Indeed, using this method arrives at 
results for the emission fraction which are very similar to 
those obtained when we use the mass deposition rate to 
define CF clusters. Th e separation is not as striking as in 
IMcCarthv et al.l (|2004T) . probably because the measured ve- 
locity dispersion can be systematically affected by substruc- 
tures. 



3.2 Magnitude Dependence 

We show in Fig.[5]the line emission strength as a function of 
K absolute magnitude for the NFPS and SDSS BCGs and 
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Table 1. Summary of Results. The first column presents the survey studied, NFPS or SDSS. SDSS (RED) includes only the red-selected 
SDSS galaxies, and SDSS (X-ray) only those with X-ray counterparts. Column 2 indicates if the results are for the BCGs or Controls. 
Column 3 shows the total number of galaxies for each sample. The number of these which are emitting is shown in column 4, and column 
5 shows the fraction of emitting galaxies. Column 6 lists the number of strongly emitting (W (Ha)> 2A) galaxies in each sample for 
which the emission line ratios are characteristic of AGN activity (usually LINER). We present the number of galaxies in our sample 
which are known to belong to a CF cluster in column 7. Column 8 shows the number of these galaxies in CF which are also emitting, 
and the final column shows the fraction of emitting galaxies in CF clusters. 



Survey Sample Total Emitting Emitting AGN or Comp. Known Emitting Emitting 

Fraction(%) W D (H«)> 2A CF CF Fraction(%) 



NFPS 


BCGs 


60 


12 


20±6 


N/A 


11 


10 


71 + 9 

"■-14 




Controls 


159 


15 


9±2 


N/A 


36 


5 


14±6 


SDSS 


BCGs 


328 


12 


13±2 


31 


N/A 


N/A 


N/A 




Controls 


526 


91 


18±2 


57 


N/A 


N/A 


N/A 


SDSS (RED) 


BCGs 


318 


35 


11±2 


25 


N/A 


N/A 


N/A 




Controls 


446 


51 


11±2 


39 


N/A 


N/A 


N/A 


SDSS (X-ray) 


BCGs 


34 


6 




3 


N/A 


N/A 


N/A 



control galaxies. The NFPS clusters are separated into CF 
and non-CF cases. In the non-CF clusters, there are no emit- 
ting BCGs, and the fraction of emission-line galaxies, shown 
in the bottom panels, is also very low for the control galax- 
ies (~ 5 %). On the other hand, in CF clusters ~ 70 % of 
the BCGs show emission, as we noticed also in the previous 
section. Moreover, the BCGs with strong emission tend to 
be the brightest galaxies, Mk < —25.5, where by definition, 
there are few control galaxies. Hence, this trend is quite dif- 
ferent from what is seen in the control sample, where almost 
all of the emission line galaxies are fainter than Mk = —25. 
For magnitudes near Mk ~ —25, where there is substan- 
tial overlap between the two populations, the emission line 
fraction of both populations are similar for non-CF clus- 
ters. The results from the SDSS, which include all clusters 
regardless of X-ray luminosity or CF status, are consistent 
with the results for the total NFPS. There are somewhat 
fewer emission-line BCGs at a given magnitude than the for 
NFPS, presumably because CF clusters make up a smaller 
proportion of the sample in the SDSS. We will say more 
about this in Section [4] 

Also highlighted on Fig.[5]and Fig.|6]are the BCGs with 
LINER-like emission. Most of the emitting BCGs are more 
characteristic of LINER emission; this is especially true of 
the less luminous BCGs. 

3.3 Dependence on Location in Cluster 

In our samples, there are some BCGs that are found close to 
the X-ray centre, while others can be found several hundred 
kpc away. In this section we will investigate whether or not 
the presence of emission lines in the BCG depends on its 
distance to the X-ray centre. 

This is illustrated in Fig. [6]where we show the line emis- 
sion as a function of cluster radius for the NFPS galaxies. All 
of the strongly emitting BCGs are within 50 kpc of the X-ray 
centre. As discussed in the previous section, these emitting 
galaxies are also usually found in a cooling flow cluster. In 
the rightmost panel we show the equivalent plot for the 34 
BCGs in the SDSS sample for which we have X-ray centres. 
Although the sample is small, our results are consistent with 
those seen in the NFPS; only those BCGs that are close to 
the X-ray centre have significant line emission. 



A Kolmogorov-Smirnov test on the H/3 distributions of 
the non-CF BCGs and the controls shows no evidence for 
a difference in the two populations. Therefore, to summa- 
rize, we conclude that an increased frequency of optical line 
emission in BCGs is observed only for those galaxies that 
lie within 50kpc of the X-ray centre of a CF cluster. For the 
remainder, the frequency of emission in BCGs is consistent 
with that observed in the control population. 

3.4 Dependence on Cluster Mass and Density 

The SDSS provides us with an optically-selected sample, 
spanning a wide range in velocity dispersion, that should 
be representative of the cluster population independent of 
X-ray properties. In this section we use this full sample (in- 
cluding those galaxies with (it' — r') < 2.2 that were ex- 
cluded when comparing directly with the NFPS) to explore 
the effect of environment on the presence of line emission in 
BCGs. 

Since the cooling flow status of a cluster might be corre- 
lated with its total mass, or central mass density, we wish to 
explore whether the trends we have observed are merely re- 
flecting a more fundamental correlation with either of these 
quantities. First, in Fig. [7] we plot the fraction of BCGs with 
Ha emission as a function of the cluster velocity dispersion. 
There is no strong trend. Rather, a fraction ~ 10-15 % of 
BCGs in each bin showing Ha emission is observed (the 
Spearman correlation coefficient is 0.17). We note that the 
frequency of radio-loud AG N is also found t o be indepen- 
dent of velocity dispersion (jBest et al.ll2006l ). The control 
galaxy population in the SDSS sample also shows no strong 
correlation between emission line fraction and group velocity 
dispersion (the Spearman correlation coefficient is 0.38). For 
the control galaxies, the overall fraction of emitting galaxies 
is somewhat higher than for BCGs, closer to 15-20%. 

Next, we examine in Fig. [7] the frequency of emis- 
sion lines in BCGs as a function of the galaxy density, as 
measured by the distance to the fifth-nearest spectroscopic 
neighbour. Even in the densest regions this corresponds to 
a smoothing scale of ~ 200 kpc, so the measurement here is 
of the total mass density on scales larger than the cooling 
radius. Generally, one would expect the central regions of 
CF clusters to be the densest environments, and if the pres- 
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Figure 5. Line Emission as a function of K-Band Magnitude. Top Panel: The equivalent width of the line emission as a function of 
K-band magnitude for the NFPS galaxies on the left, and for the SDSS galaxies on the right. The NFPS galaxies are in subsets: a) All 
galaxies, including those in clusters where the CF status is unknown, b) those which are in known CF clusters, and c) those which are in 
known non-CF clusters. Filled circles represent BCGs, Xs for control galaxies, and open circles indicate LINER emission (for the SDSS 
sample only). The dashed line represents the cut between emitting and non-emitting galaxies. The BCGs for Abell 780 and Abell 1795 
have unreliable H/3 equivalent widths, but nonetheless strong emission, and are represented by lower limits at W > 3.3 A and > 3.7 A, 
respectively. Bottom Panel : Using the same subsamples as the top panel, we plot the fraction of line emitting galaxies as a function of 
K-band magnitude. The solid line represents the BCGs, and the dotted line the control galaxies. 
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Figure 6. Line Emission as a function of distance from the cluster X-ray centre. Subsample definitions and symbols are the same as in 
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Figure 7. Left Panel: The fraction of Ho emitting galaxies in the SDSS as a function of group velocity dispersion. The solid line refers to 
BCGs and the dotted line are the control galaxies. This adaptive histogram contains 85 galaxies per bin for the BCGs and 100 galaxies 
per bin for the controls. Right Panel: The fraction of Ha emitting galaxies in the SDSS as a function of the galaxy density. This adaptive 
histogram contains 85 galaxies per bin for the BCGs and 95 galaxies per bin for the controls. 



ence of a cooling-flow was correlated with the mass or galaxy 
density on larger scales our previous results would lead us to 
expect the most star formation occurring in the most dense 
regions. On the contrary, we observe a clear correlation with 
the emission line fraction decreasing with increasing number 
density, and so BCGs are less likely to show emission lines 
if they are found in the densest regions of clusters (in both 
cases of the BCG and the controls, the Spearman test yields 
a correlation coefficient of ~ 0.9, and the correlation is sig- 
nificant at the ~ 95% confidence level). At all densities, the 
fraction of emission-line galaxies is higher in the control pop- 
ulation than the BCG population. The fact that we observe 
enhanced emission for those galaxies within the smaller 50 
kpc scale of the X-ray centre of CF clusters (Section 13. 3p 
is therefore very likely related directly to the presence of 
cooling, X-ray gas on small scales, rather than the overall 
gravitational potential. 



4 DISCUSSION 

The overall fraction of BCGs with emission lines is ~ 13 % 
in the SDSS, an d ~ 20 % in the NFPS. The latter is in good 
agreement with ICrawford et al.l l|l999l ) who find a fraction 
of 27 % emission-line BCGs in their sample of X-ray se- 
lected clusters. We do not know what the CF fraction is 
in an optically-selected sample, or as a function of mass. 
Nonetheless, since the fraction of massive cluste rs that host 
a cooling flow is likely no more than about 50 % l|Peres et al.l 
ll99Sl ; lMcCarthv et al.ll200l ; IChen et alj|2007h . and CF clus- 
ters are systematically overluminous, we attribute the factor 
of two difference in emission fracti on between our two s ur- 
veys to the fact that NFPS (and ICrawford et al.l \l99$ ) is 
X-ray selected, and therefore biased toward CF clusters. 

There is currently much observational and theoret- 
ical work exploring the possible feedback between a 



central galaxy's AG N, current s tar fo r mation, and the 
cooling X-ray gas (ISilk fc Reesl 1 19981 : iBest et all 120061; 
I Croton et ail l200f]|; iBower et all 120061 : ISiiacki fc Springell 
boodlDe Lucia fc Blaizotll2006l . e.g.). AGN in the BCGs are 
thought to play an important role in suppressing cooling, 
and lowering mass deposition rates (and hence star forma- 
tion rates). Our main result is that the majority of line emit- 
ting BCGs are positioned close to the X-ray centre of clusters 
classified as hosting a cooling flow directly from the mass de- 
position rates (as defined in Section[2TT} . Furthermore, this 
result holds when we classify cooling flow clusters instead 
by their excess X-ray luminosity relative to other clusters 
with similar dynamical mass (a definition we explored at 
the end of Section l3~Tj) . Similar conclusions, based on smaller 
samp l es, have been reached by other s (e.g. Johnstone et all 
19871: iMcNamara fc Q'Connelllll989l : ICrawford et all Il999l : 



Raffertv et al.ll2006t ). Importantly, we have shown that, in 



the absence of a cooling flow, emission lines are rare in 
BCGs, regardless of the mass density or velocity dispersion 
of the cluster. Moreover, a control population of similarly 
bright, centrally located cluster galaxies does not exhibit 
an increased frequency of emission lines in CF clusters. We 
therefore conclude that the observed emission (arising from 
star formation and/or an AGN) is directly related to the 
presence of cooling gas. 

For the most part, we have not concerned ourselves with 
the origin of the observed emission, since starburst galax- 
ies and optically-selected AGNs are probably closely linked 
iKauffmann et all |2003). However, it is wo rth investigat- 
ing this further here. Crawford et all (|l999h find that the 
very strongest H a emitters have star formatio n-like emission 
line ratios, while Ivon der Linden et all l|2006l ) find that most 
BCGs, which display weaker Ha emission, are more charac- 
teristic of LINER activity. The six NFPS BCGs that have 
[Nil] /Ha line ratios available from ICrawford et all (|l999l ) 
lie in the regime straddling LINERS and Seyferts. Of the 
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six emitting BCGs in the SDSS-X-Ray sample, the three 
stronger Ha emitters are non-AGN, which means that they 
are likely composite, whereas the somewhat weaker emitters 
are classified as AGN. Thus the emission-line BCGS found 
near the X-ray centre of galaxy clusters appear to be a het- 
erogeneous class of objects. 

In the unrestricted SDSS sample, we find only ~ 13 % of 
BCGs have line emission, compared t o ~ 18 % for controls. 
This e xhibits the same trend as in the lvon der Linden et all 
(2006) study of SDSS C4 cluster galaxies, where detectable 
emission line luminosity measurements (with S/N >3) are 
found in 30 % of BCGs, as opposed to 40 % for the other 
massive, central, galaxies (ie. controls). Our numbers can- 
not be compared directly, as their cut includes more weakly 
emitting systems as emitters than does our more strict cri- 
terion of Wo > 2 A. If we relax our definition of an emitting 
galaxy being one with W > A, we find the percentage 
of line emitters increases to 27 % for BCGs, and 38 % for 
controls. 

Our Fig. [7] reiterates the above point, that in the SDSS 
sample, the BCGs have a lower fraction of line emitting 
galaxies than the controls. This is also consistent with what 
iBest et all (|2006h find: that emission line AGN activity is 
suppressed in the galaxies near the centre of the cluster with 
respect to the other massive cluster galaxies. They point out 
that emission line AGN and radio-loud AGN are partly in- 
dependent populations, but that BCGs with emission-line 
AGN are more likely to host radio-loud AGN than other 
galaxies. Best et al. also find that radio-loud AGN are pref- 
erentially found in BCGs within 0.2 r2oo of the centre of 
the cluster and that there is a strong dependence on galaxy 
stellar mass, with ~ 40 % of the most massive galaxies show- 
ing radio loud AGN emission. Therefore it is plausible that 
our emission-line BCGs in the NFPS (which are found at the 
centre of CF clusters) and in the SDSS (many of which show 
LINER-like activity), host radio-loud AGN. We therefore 
matched our NFPS BCGs to radio sources from the Faint 
Images of the Ra dio Sky at Twenty-Centimeters survey 
JWhite et al.lll997l. hereafte r FIRST) and NRAO VLA Sky 
Survey l|Condon et al.|[l998l. hereafter NVSS). We find that 
indeed, all 12 of our H/3 emitting BCGs are radio sources. 
Furthermore, for the 50 BCGs in the FOV of the surveys, 
28 of our BCGs have radio counterparts (~ 56±11 %). 
Thus ~ 42±11 % of non-H/3 emitting BCGs have radio 
emission. Again, the cooling flow status of the cluster ap- 
pears important: ~ 71+^ % of CF BCGs are radio sources, 
and ~ 80t^ % of central CF BCGs are radio sources. For 
non- CF BCGs only ~ 32±^ % have radio emission. For the 
controls, 36 of 144 galaxies in the FOV of the surveys are 
radio sources (~ 25±4 %). 

Recently, iLin fc Mohrl (|2006^ have studied the radio- 
loud properties of a large X-ray selected cluster sample, and 
find the overall radio-loud fraction to be ~ 35 % for PC- 
band selected BCGs within r2oo of the cluster center (com- 
pared with ~ 20 % for bright cluster galaxies, excluding 
the BCGs). This is in reasonable agreement with our re- 
sults, although the radio-loud BCG fr action we measure 
is somewhat higher. ILin fc Mohrl | 2006) also find a strong 
trend with cluster mass (inferred from the X-ray luminos- 
ity), though the strength of the trend is sensitive to the radio 
power limit and K-band luminosity of the galaxy. Our NFPS 
sample is too small to robustly identify such trends, but we 



note that the fraction of BCGs with radio sources that are in 
clusters with a velocity dispersion < 600km/s is 54±15 %, 
quite similar to that for BCGs in clusters with a velocity 
dispersion > 600km/s, 58±15 %. In the optically-selected 
SDSS, we find no significant trend with velocity dispersion, 
but an overall fraction of ~ 40 % of the BCGs with radio 
emission, and ~ 25 % of the controls. Because of the differ- 
ence in sample selection, and the cluster mass estimators, we 
d o not conside r this a serious discrepancy with the results 
of ILin fc Mohrl (|2006l ). 



5 CONCLUSIONS 

We have used two large, homogeneous galaxy cluster surveys 
to investigate the incidence of optical emission lines amongst 
BCGs. The NFPS consists of 60 BCGs in X-ray selected 
clusters, while the SDSS sample is a larger, optically-selected 
sample of 328 BCGs. From these data, we are able to draw 
the following conclusions: 

• Of the 10 BCGs that lie within 50kpc of the peak of 
X-ray emission in a cluster with evidence for a significant 
cooling flow, all show optical emission lines. Moreover, of the 
12 BCGs that show emission, all are located within 50kpc of 
the X-ray emission peak, all are radio sources, and none are 
in clusters of known non-CF status (10 are in CF clusters, 
and 2 are in clusters with unknown CF status). 

• Excluding the special circumstances noted above, the 
fraction of BCGs that exhibit optical emission lines is ~ 10- 
20 %, and is always comparable to or lower than the fraction 
for control galaxies with a similar luminosity and environ- 
ment. 

• For optically selected cluster samples, which are domi- 
nated by non-CF clusters, the fraction of BCGs with emis- 
sion does not correlate strongly with cluster mass or galaxy 
density. 

We have therefore demonstrated a direct connection 
between the presence of cooling gas, and enhanced optical 
emission in a centrally located galaxy. It would be very useful 
to obtain pointed X-ray observations of those SDSS clusters 
in which we have found a BCG with Ha emission, to deter- 
mine if this correlation holds in optically-selected samples. 
These clusters would also be potentially interesting for ob- 
servation with Chandra to observe the X-ray emission mor- 
phology, as other massive clusters with Ha emitting BCGs, 
such as Abell 426 and Abell 1795 show remarkable X-ray 
morphology, such as X-ray holes, and cooling tails. 
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Table 2. Tabic of NFPS BCGs. The cluster name is shown in column 1. Column 2 and 3 give the position of the cluster X-ray centre. 
The cluster redshift is given in column 4, the cluster velocity dispersion and r200 are given in columns 5 and 6. Column 7 gives the cooling 
flow mass deposition rate (MDR) in (Mfl/yr). Colum n 8 gi ves the cooling flow stat us of the cluster. Colum n 9 gi ves the reference fo r 
the MDR/CF status: b s tands for lBirzan et al.l (1200411 f for l Fuiita et al.l <2006ft , i for I Johnstone etal] j2005T). h for lHenrv et al.l <2004l) 
k forlKempner &c Davidl d2004h m for lMcCarthy et all (I2004T) , p for lPeres et all l ll99Sfi . r for lSharma et al.l \2004 ), s for lSanderson et al.l 



(20061. v for lKanov etaiT l2006'). w for I White! <200C|) , and o for other. Column 10 is the X-ray luminosity in units of 10 44 erg/s, column 
11 gives the name of the BCG, column 12 is the BCG K-band magnitude, column 13 gives the distance between the BCG and the cluster 
X-ray centre. Columns 14 gives the H/3 equivalent width, and column 15 gives the error, which is the noise in the line-free regions of the 
absorption-corrected spectrum. 
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